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Abstract. Kirchhoff 's matrix-tree theorem states that the number of spanning trees of a graph 
G is equal to the value of the determinant of the reduced Laplacian of G. We outline an efficient 
bijective proof of this theorem, by studying a canonical finite abelian group attached to G whose 
order is equal to the value of same matrix determinant. More specifically, we show how one can 
efficiently compute a bijection between the group elements and the spanning trees of the graph. The 
main ingredient for computing the bijection is an efficient algorithm for finding the unique G-parking 
function (reduced divisor) in a linear equivalence class defined by a chip-firing game. We also give 
applications, including a new and completely algebraic algorithm for generating random spanning 
trees. Other applications include algorithms related to chip-firing games and sandpile group law, as 
well as certain algorithmic problems about the Riemann-Roch theory on graphs. 



1. Introduction 

1.1. Overview. Every graph G has a canonical finite abelian group attached to it. This group has 
appeared in the literature under many different names; in theoretical physics it was first introduced 
as the "abelian sandpile group" or "abelian avalanche group" in the context of self-organized critical 
phenomena ( [H EOJ [2T] ) . In arithmetic geometry, this group appeared as the "group of components" 
in the study of degenerating algebraic curves f[3Uj). In algebraic graph theory this group appeared 
under the name " Jacobian group" or "Picard group" in the study of flows and cuts in graphs ( [3] ) . 
The study of a certain chip-firing game on graphs led to the definition of this group under the name 
"critical group" (0Q2]). 

The order of this group is equal to the value of the determinant of the reduced Laplacian of G 
(see, e.g. Lemma [2 .3p . We know from Kirchhoff 's famous matrix-tree theorem that the value of the 
same determinant gives the number of spanning trees of the graph (|27J). So, one might wonder 
whether there is a nice and explicit bijection between the elements of the group and the spanning 
trees of the graph; existence of such a bijection would independently prove the matrix-tree theorem, 
and might have other algorithmic consequences. It is rather clear that such a bijection cannot be 
fully canonical, as that would imply a particular spanning tree is distinguished, and corresponds to 
the identity of the group. Therefore, one needs to make some choices to be able to write down a 
bijection. If one fixes a vertex q, then there is a canonical representative for each element of the 
group, called the G-parking function (based at q) or q-reduced divisor (see, e.g., |21l 1191 [3"Tl [5] , or 
Proposition 12 . 5[) . The main result of this paper is an algorithm to find this canonical representative 
efficiently. Once we have this, we can use one of the bijections in [12\ \TE[ [TBI E] to find the 
corresponding spanning tree. 

This easy-to-compute bijection from the "Jacobian group" to the set of spanning trees of the graph 
gives new insight into the matrix-tree theorem. It also may lead to new algorithms for graphs and 
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their spanning trees. For example, as an immediate corollary, if one picks a random element of the 
group - which is a trivial task - and constructs the bijection, then one gets a random spanning tree. 
This yields a new and completely algebraic approach for sampling a random spanning tree from the 
set of spanning trees of the graph. With this approach, it is very easy to sample multiple spanning 
trees with certain "joint" distribution. We believe there may be other algorithmic applications for 
this bijection. 

The theory of G-parking functions relates nicely to a number of problems in theory of chip-firing 
games , and to the Riemann-Roch theory on finite graphs. As a result of our main algorithm, we 
are able to give efficient algorithms for a number of problems in these areas. The question of finding 
reduced divisors (G-parking functions) was first posed by Henrik Lenstra to Baker and Norine, 
in connection with Riemann-Roch theory on finite graphs in [5], through private communication. 
The problem of finding the "sandpile prediction" in polynomial-time was listed as an open problem 
by Laszlo Babai in [2] (§10.3). Our main algorithm also settles this open problem, as the sandpile 
predictions and critical configurations are, in a sense, dual to the reduced divisors; see Remark l3.141 

1.2. Some related work. The term "G-parking functions" was first introduced in [31] . The 
reason for this terminology is that they can be considered as a natural generalization of "parking 
functions" . The theory of parking functions was first considered in connection with hash functions 
(|28j). The original problem was phrased in terms of cars and parking spots. The theory of parking 
functions has since been developed, with links many different areas including priority queues (|22j), 
representation theory ([24]). and noncrossing partitions (|32j). Although the explicit definition was 
first given in [31], the concept had appeared (sometimes in disguise) in many previous work (see, 
e.g., [201 ED El ED IS] ). The fact that G-parking functions provide a canonical representative for 
each element of the Jacobian appear implicitly in [2D US], and explicitly in [5T| l5j. 

The problem of giving explicit bijection between G-parking functions and spanning trees of graph 
is studied extensively in the literature (see, e.g., [121 HHJ [I6j [7] ) . 

The relationship between chip-firing games and the Jacobian group is studied in [SI [10] . Some 
algorithmic aspects of the chip-firing games are studied in [311 [HI [33] ; see §A.2l for a discussion on 
how [34] relates to our work. 

The Uniform Spanning Tree (UST) problem has been extensively studied in the literature and 
there are two known types of algorithms: determinant based algorithms (e.g. [231 fl7l I29j). and 
random walk based algorithms (e.g. [151 [p [35]). 

More related work is given in the paper. 

The paper is structured as follows. In £j2]we provide the relevant definitions and tools. In £j3]we 
first state and prove Dhar's algorithm which can efficiently check whether a given divisor is reduced. 
Then we give our main algorithm for computing the reduced divisors and prove its correctness and 
efficiency. We then give the explicit and efficient bijection between the Jacobian group and the set 
of spanning trees of the graph. In $A]we outline some applications of our main results, including a 
new approach to the Uniform, Spanning Tree problem, efficient computation of the group law in the 
critical and sandpile group, efficient algorithm for the chip-firing game of Baker and Norine, and 
the relation to the Riemann-Roch theory on finite graphs. 

2. Definitions and background 

2.1. Notation and Terminology. Throughout this paper, a graph means a finite, unweighted 
multigraph with no loops. All graphs are assumed to be connected. For a graph G, the set of 
vertices is denoted by V(G), and the set of edges is denoted by E(G). Throughout this paper, n 
and m denote the number of vertices and edges, respectively. 
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Let Div(G) be the free abelian group generated by V(G). One can think of elements of Div(G) 
as formal integer linear combination of vertices 

Div(G) = { a v(v) ■■ a v e Z} . 

By analogy with the algebraic curve case, elements Div(G) are called Divisors on G. For a divisor 
D, the coefficient a v of (v) in D is denoted by D(y). 

Ai(G) = Hom(V(G),Z) is, by definition, the abelian group consisting of all integer- valued func- 
tions on the vertices. For A C V(G), xa £ M(G) denotes the {0, l}-valued characteristic function 
of A. Note that {X{v}}v&v(G) generates M.(G) as a Z-module. 

The Laplacian operator A : M.{G) — > Div(G) is given by the formula 

A(/)= A v(f)(v), 
ueV(G) 

where 

A„(/) = Y " /W) • 

{v,w}£E(G) 

Let {vi,... ,v n } be an ordering of V(G). With respect to this ordering, the Laplacian matrix 
Q associated to G is the n x n matrix Q = (qij), where qu is the degree of vertex V{, and — qij 
(i 7^ j) is the number of edges connecting and Vj. It is well-known (and easy to verify) that Q is 
symmetric, has rank n — 1, and that the kernel of Q is spanned by 1, the all-one vectorlj (see, e.g., 

BUM). 

Using the ordering {v\ , . . . , v n } of vertices, there exist isomorphisms between abelian groups 
Div(Gr), A4(G), and the Z-module of n X 1 column vectors having integer coordinates. Under 
these isomorphisms the Laplacian operator A : Ai(G) — > Div(G) coincides with the Z-module 
homomorphism Q : Z n — > Z n . More specifically, if [D] denotes the column vector corresponding to 
D G Div(G), and [/] denotes the column vector corresponding to / € A4(G), then [A(/)] = Q[f]. 

2.2. Chip-firing games on graphs. Following [5j we define an equivalence relation ~ on the 
group Div(G) as follows. 

Definition. For D\,D% G Div(G), D\ ~ L> 2 if and only if Z?i - D 2 is in the image of A : M(G) — > 
Div(G). 

It is an easy exercise to show that this is indeed an equivalence relation. This equivalence 
relation is closely related to notion of "chip-firing games" played on the vertices of the graph (see, 
e.g., [13 [9j Uni E] ) . For a given divisor D € Div(G), the integer D{y) can be viewed as the amount 
of dollars that is assigned to the vertex v, and D can be viewed as a configuration of the economy 
in the society modeled by the graph. If D(v) < then the vertex v is said to be in deb^. A move 
in this "one-player game" consists of choosing a vertex and having it borrow one dollar from each 
of its neighbors, or give ("fire") one dollar to each of its neighbors. 

The following easy lemma is proved in [5]. 

Lemma 2.1. For Di,D 2 £ Div(G), D\ ~ D 2 if and only if starting from the configuration D\ one 
can reach to the configuration D 2 , through a sequence of moves. 



We again emphasize that G has no loops. 
2 Since negative assignments to vertices are permitted, we prefer to use the term "dollar" instead of "chip". This 
term was first used by Biggs in [9l 110]. 
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2.3. The Jacobian of a finite graph. Consider the group homomorphism deg : Div(G) — > 
Z defined by deg(-D) = X^eV(G) D(v). Denote by Div°(G) the kernel of this homomorphism, 
consisting of divisors of degree zero. Denote by Prin(G) the image of the Laplacian operator A : 
M{G) -> Div(G). 

Lemma 2.2. Prin(G) C Div (G), and both Prin(G ! ) and Div°(G) are free 7L-modules of rank n — 1. 

Proof is easy and is left as an exercise. 

As a corollary, the quotient group 

Jac(G) = Div°(G)/Prin(G) 

is well-defined and is a finite abelian group. Following [3], it is called the Jacobian of G. By the 
fundamental theorem of finitely generated abelian groups, this group is isomorphic to direct sum of 
finite cyclic groups. 

If we choose an ordering of vertices it is easy to see that Jac(G) is the torsion par10 of the cokernel 
of the Z-module homomorphism Q : Z n — ► Z n . The following result is an easy consequence of this 
observation. 

Lemma 2.3. The order o/Jac(G) is equal to the value of the determinant of the reduced Laplacian 
of G. Invariant factors and generators of Jac(G) can efficiently be computed by finding the Smith 
normal form of the Laplacian matrix of the graph G. 

Remark 2.4. Of course the value of the determinant of the reduced Laplacian of G gives the number 
of spanning trees of G. This is Kirchhoff's matrix-Tree Theorem (|27j). We will not assume the 
knowledge of this famous theorem and, instead, we will use Lemma [2.3l to give a new bijective proof, 
with computable bijection, of the matrix-tree theorem in £|3.5l Indeed, one can think of Lemma [2.3I 
as the generalization of Kirchhoff's theorem; the matrix-tree theorem only computes the order of 
the group and forgets the structure!. 

2.4. Reduced divisors or G-parking functions. For A C V(G) and v G A, let outdeg^(u) 
denote the number of edges of G having v as one endpoint and whose other endpoint lies in V(G)\A. 

Definition. Fix a vertex q G V(G). A function / : V(G)\{q} — ► Z is called a (/-reduced function 
(or G-parking function based at q) , if it satisfies the following two conditions: 

• f(v) > for all v G V(G)\{q}. 

• For every non-empty set A C V(G)\{q}, there exists a vertex v G A such that f(v) < 
outdeg A (u). 

This definition might seem unmotivated and rather technical for the moment. Its significance 
will only be clear after Proposition I2.51 

In this paper, following [5], we use the term q-reduced. Moreover, a divisor D G Div(G) is called 
q-reduced if the map v i— ► D(v), defined for v G V(G)\{q}, is a G-parking function (relative to the 
base vertex q). 

The significance of (/-reduced divisors is mainly due to the following proposition. This was 
discovered by several different authors (see, e.g., (2TJ [T21 El], [5] ) 

Proposition 2.5. Fix a base vertex q G V{G). Then for every D G Div(G), there exists a unique 
q-reduced divisor D' G Div(G) such that D' ~ D. 

Remark 2.6. (i) None of the proofs in the literature gives to an efficient algorithm for computing 
such D' ~ D. This work is done in the current paper. 



'The full cokernel group gives the full "Picard group" Pic(G) = Div(G)/Prin(G), which is isomorphic to Z©Jac(G) 
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(ii) There is a close relationship between g-reduced divisors and g-critical configurations (or 
sandpile predictions); see Remark 13.141 for more detail. 

(iii) It follows from Proposition 12.51 that the Jacobian of G induces a group structure on the set 
of q-reduced divisor, or on the set of g-critical configurations (see £ jA.2l for more about this 
observation). However, the group law on these sets will be somehow "artificial". 

2.5. Generalized inverses. A matrix can have an inverse only if it is square and its columns (or 
rows) are linearly independent. But one can still get "partial inverse" of any matrix. 

Definition. Let A be a matrix (not necessarily square). Any matrix L satisfying ALA = A is 
called a generalized inverse of A. 

It is somehow surprising that for every matrix A there exists at least one generalized inverse. In 
fact, more is true; any matrix has a unique Moore-Penrose pseudoinvers^. 

Let Q be the Laplacian matrix of a connected graph. Since it has rank n — 1 it cannot have an 
inverse. But there are many ways to obtain generalized inverses. Two examples are given below: 

Example 2.7. Fix an integer 1 < % < n. Let Qi be the (rt — 1) x (n — 1) matrix obtained from Q 
by deleting i th row and i th column^. Then Qi is a full rank matrix and has an inverse Q~ l ■ Let Lu\ 
be the n x n matrix obtained from Q^ 1 by inserting a zero row after (i — l) th row and inserting a 
zero column after (i — l) th column. Then is a generalized inverse of Q. One can check 

QL(i) = I + R(i) 

where / is the identity matrix, and has —1 entries in i th row and is zero everywhere else. As 
R(i)Q = 0, we get QL^Q = Q. 

Example 2.8. Let J be the n x n all one matrix. Then Q + -J is nonsingular and Q + = 
(Q + ^J) 1 — \ J is a generalized inverse of Q. In fact it is the unique Moore-Penrose pseudoinverse 
of Q. It is easy to check QQ + = Q + Q = I — ^ J. 

3. Algorithm for computing reduced divisors 

3.1. Dhar's algorithm: the decision problem. Given a divisor on the graph G, definition 
suggests that the verifier needs to check whether there exists a vertex v £ A such that D(v) < 
outdeg^f) for all subsets A C V(G)\{q}. But, in fact, there is an elegant algorithm to do this 
quickly, which is called the Dhar's algorithm (after Dhar [20J). Both the algorithm and its proof of 
correctness are given here. 

Algorithm. (Dhar's algorithm) 

Input: A divisor D € Div(G), and a vertex q £ V{G). 

Output: TRUE if D is ^-reduced, and FALSE if D is not ^-reduced. 

Let = V(G) and vq = q. 

(0) Check D(v) > for all v G V(G)\{q}. 

(1) For 1 < i < n — 1, let A^ = Ai^i \ Find a V{ & A{ with D(vi) < outdeg^. {v{) - If no 
such Vi exists, output FALSE and stop. 

(2) Output TRUE. 



Moore-Penrose pseudoinverse of A is a generalized inverse of A with three extra properties; see [6] for an extensive 
study of the subject. 

^Qi is called a reduced laplacian of the graph. 
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Analysis of the Dhar's algorithm. If the algorithm outputs FALSE, D cannot be g-reduced by 
definition. Now assume the algorithm outputs TRUE. We need to show that for an arbitrary 
subset A C V(G)\{q}, there exists a vertex v £ A such that D(v) < outdeg^(v). Let vp. be the 
minimal element of A, with respect to the ordering induced by the algorithm. Then A C Aj. and 
therefore outdeg^ (vk) < outdeg j4 (vfe). But we know D(vk) < outdeg^ fe (?;fc) by construction. Hence 
D(vk) < outdeg yl (ufc) and Vk is the desired element in A. 

The Dhar's algorithm clearly runs in time 0(n 2 ) . □ 

3.2. The main algorithm: the search problem. Fix a base vertex q £ V{G). Here is the 
problem we want to solve. 

Given a divisor D E Div(G), we want to find the unique (/-reduced divisor D' ~ D. 

A basic idea is to try to make the Dhar's algorithm to output TRUE; anytime it is about to 
output FALSE, "modify" the divisor to avoid it. However, there is a major problem with this naive 
idea; this process can take exponential time in the size of input. Of course, one also needs to figure 
out how to "modify" the divisor when the Dhar's algorithm is about to fail. 

For us, the multigraph can be described by the number of vertices n and edges m in the graph. 
A divisor D can be presented by roughly YlveV(G) ^og(D(v)) bits, which is easily seen to be less 
than n ■ log(deg(-D)) bits. 

Remark 3.1. In fact, it is not too hard to find an algorithm that finds the g-reduced divisors in time 
that runs in polynomial in n, m, and deg(D) using the existing techniques appeared in [T3l l33j 125]. 
But this is exponential in the "size of input" ! Our algorithm will have a running time polynomial 
in m and n, and there will be some log(deg(-D))-bit computations involved. 

We will now give an algorithm for finding the reduced divisor and prove its correctness in ^3.31 In 
§3.41 we give a deterministic upper bounds for its running time which show the algorithm is indeed 
efficient. 

Notation. |_"J denotes the floor function. For a vector X, we denote by [X\ a vector whose 
entries are the floor of entries of X. Recall from §2.11 that xa £ M(G) denotes the {0, l}-valued 
characteristic function of A C V{G). To simplify the presentation, we pick an ordering of vertices, 
and use [•] for column vector presentation of divisors and functions. 

Algorithm. (Finding the reduced divisor) 

Input: The Laplacian matrix Q of graph G. A divisor D G Div(G), and a vertex q £ V(G). 
Output: The unique g-reduced divisor D' ~ D. 

Step 1. Find the generalized inverse L = of Q, as in Example 12.71 Compute the divisor 
[D'] = [D]-Q[L {q) [D}\. 

Step 2. Find a vertex v / q with D'(v) < and substitute [D'\ with [D'\ + Q[x{v}] (i- e - have v 
borrow one dollar from each of its neighbors). Repeat Step 2 until no such "negative" vertex 
v ^ q exists anymore. 

Step 3. Now let A = V(G) and v = q. 

(1) Set i = 1. 

(2) While i < n - 1, let A { = A^ \ 

* If there exists a vi G Ai such that D'{vi) < outdeg^. (vi): increase i by 1 and go 
back to (2) 

* Else, compute 

k= mm [D'(vi)/out&eg A (vi)\ 

outdeg A .(t)i)^0 
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and replace [D 1 ] with [D'\ — Q[k ■ xaJ (i- e - have all vertices in Ai fire enough 
times to get D'(v{) < outdeg^. (vi) for at least one vertex in A{). Then go back 
to (1). 

3.3. Correctness of the algorithm. Here is the high-level idea of the the steps in the algorithm 

• The task of Step 1 is to bound the values of the divisor on vertices v ^ q. Bounding these 
values is important for bounding the running time of Step 2 and Step 3. 

• The task of Step 2 is to make the values of the divisor on v ^ q non-negative while still 
keeping the values bounded. 

• The task of Step 3 is to "force" the divisor to pass the Dhar's algorithm. 

Assume for the moment that the algorithm actually terminates and produces an output. It is 
easy to see that the output would be equivalent to D. Also, the output would certainly pass the 
Dhar's algorithm and therefore is g-reduced. Therefore, for the correctness of the algorithm, we 
only need to show that it always terminates and produces an output. 

Step 1. The following proposition clarifies the task of Step 1. 

Proposition 3.2. // [D 1 ] = [D] - Q[L (q) [D}\, then \D'(v)\ < deg{v) for all v / q. D'(q) is such 
that deg(D') = deg(L>). 

Proof. Recall from Example 12.71 that QL^ = I + Rtq), where I is the identity matrix and has 
— 1 entries in g th row and is zero everywhere else. Therefore [D] = QL^[D] + deg(D) • e q , where 
e q is the column vector which is zero everywhere except at position q, which is 1. Now 

[D'] = [D]-Q[L {q) [D}\ 

= Q(L (q) [D]-[L (q) [D}\)+deg(D)-e q 
= Qf + deg(D) ■ e q 

where f = L^[D] — [L^[D]\ is a vector with entries from the interval [0,1). It is now easy to 
show that the absolute values of the entries of Qf are bounded by the degree of the corresponding 
vertices. The second statement is trivial. □ 

Step 2. Perhaps it is not even clear that Step 2 will ever terminate, let alone in a polynomially 
bounded number of iterations. The following Proposition justifies this fact. 

Proposition 3.3. Given a divisor and a fixed vertex q, there exists a sequence of borrowing by 
vertices in debt (i.e. having negative numbers) which takes all vertices v ^ q out of debt (i.e. there 
exists a sequence of vertices such that Step 2 of the algorithm terminates). 

Moreover, any sequence of borrowing by vertices with in debt terminates in the same number of 
steps, with the same terminal configuration. 

Proof. This is follows from the proof of Lemma 5.3 and discussion in page 24 of [5]. 

□ 

Remark 3.4. Step 2 has very nice features; by Proposition 13.21 < deg(w) for all v ^ q. So 

for any v ^ q with D'{v) < only one borrowing is needed to make the vertex positive. Moreover, 
the resulting positive number will clearly be less than the vertex degree again. This fact, together 
with Proposition 13.21 guarantees that the output of Step 2 satisfies < D'(v) < deg(t> ) for all v ^ q. 

Step 3. In Step 3, the algorithm forces the divisor to pass the Dhar's algorithm. If for a subset 
Ai no Vi G Ai with D'(vi) < outdeg^. (vi) exists (i.e. a counter-example for being (/-reduced is 
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found), the algorithm computes 

(3.5) k = min [D' (vi)/outdeg A (vi)\ . 

outdeg A . (vj)^O 

Since 

k < D'(vi)/outdeg Ai (vi) 

or 

< D (vi) — k ■ outdeg^. (fj) 

for all V{ £ Ai with outdeg^. (t>j) ^ 0, replacing [D'\ with [D'\ — Q[k ■ XA t ] wm n °t make any vertex 
negative. Let w £ Ai be the vertex that achieves the equality f)3.5[) . Then 

D'(w)/outdeg A .(w) - 1 < k 

or 

D'(w) — k • outdeg^. (w) < outdeg^. (w) . 

This means that replacing [D'\ with [D'\ — Q[k-XA l ] will result in D'(w) < outdeg^. (w). Algorithm 
then restarts the Dhar algorithm. 

It is easy to see Step 3 eventually terminates; as vertex q never "fires" , it should eventually stop 
receiving from its neighbors at some point. This, in turn, means that its neighbors will stop firing 
at some point and therefore, eventually, will stop receiving from their neighbors. Iterating this 
argument shows that the whole process will stop at a finite number of steps. 

3.4. Bounds on the total number of chip-firing moves. It seems to be a hard problem to find 
a good estimate of the running time of the algorithm, as it depends on the structure of the graph 
and the given divisor in a rather complicated way. In practice, the algorithm seems to run much 
faster than what is predicted by the upper bound given in this section. One reason for this is that 
the technique is rather blind to the algorithm itself. This point will be clear shortly. 

The following lemma is a generalization of the standard Cauchy-Schwarz inequality which can be 
proved similarly. 

Lemma 3.6. For any positive semidefinite matrix A with largest eigenvalue n, and for all vectors 
x and y, 

\x T Ay\ < 7/[|ac[| 2 [|y[|2 

Let A = L q as in Example 12.71 Non-zero eigenvalues of L q are the reciprocal of eigenvalues of 
the reduced Laplacian matrix Q q . So is positive semidefinite and its the largest eigenvalue is 
I/A2, where A2 is the smallest eigenvalue of Q q . 

Let y = Qx. Multiplying both sides by L g , we get L g y = (I + i?^)x, or 

(3.7) L g y = x-x(g)l 

where 1 denotes the all-one vector. Then l T L q y = X]i x W — 71 ' x (?)- Now using Lemma [ 



^x(i) - n ■ x(g)| < ^-||l|| 2 ||y||2 = ^IMb 



A 2 A 2 

t 

In particular, if x(g) = we have 



(3.8) lE x «I^H2 



In fact something stronger is true. 
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Proposition 3.9. Let y = Qx where Q is the Laplacian matrix of the graph G. If x(q) = then 



, „y|li 



where \\y\\' 2 = ySi^ g yW 2 (i- e - one can assume y(q) = when computing \\y\\2 in inequality 



Proof. This follows from the above discussion, together with the following observation; in (|3.7p . 
since the g th column of is zero, one can freely change the value of y(q). Choosing y(q) = we 
are done. 

□ 



If D\ ~ D2 via a sequence of moves that q never participates, then there exists a vector x 
with x(q) = 0, such that [D{\ — [D2] = Qx. Moreover, if all moves are borrowing or all moves 
are firings, then all x(i)'s all the same sign, and | X]i x WI precisely counts the number of vertex 
moves. By Proposition 13.91 an upper bound on the number of vertex moves is given by | ^i x (^)l — 
Apll[-^l] ~~ [-^2] II 2- Since || • H2 < || • ||i for (finite-dimensional) vectors we have the following. 

Proposition 3.10. Assume D2 is obtained from D\ with a sequence of moves so that 

• q never makes a move. 

• all moves are of the same type; vertices only fire or only borrow throughout the process. 
Then an upper bound on the number of vertex moves is given by 

^(ii^iM + iim) 
A2 

where ||y||^ = £^Jy(*)l- 

Corollary 3.11. Each of Step 2 and Step 3 in the algorithm terminates in at most vertex 
moves. 

Proof. The input and output of Step 2 satisfy |-D'(f)| < deg(u) (for v 7^ q) by Proposition 13.21 and 
Remark 13.41 Therefore || • \\\ of both input and output is less than Ylv^g I deg(t>)| < 2m. Since 
vertices only borrow in Step 2 and q never makes a move, Proposition 13.101 applies. Step 3 follows 
similarly because the input and output of Step 3 satisfy |-D'(i>)| < deg(u) (for v / q) by Remark 13.41 
and definition of (/-reduced divisor. Vertices only fire in Step 3 and q never makes a move. □ 

Remark 3.12. (i) Step 3 in the algorithm performs many moves at the same time (e.g., set 
firings can be considered as combination vertex firings). So the actual number of iterations 
in Step 3 is, in practice, a fraction of what is given in Corollary 13.111 Similarly, it is possible 
to modify Step 2 to have combination moves. We omit the details here, 
(ii) In Step 1, computing the generalized inverse is the expensive part, which takes time 
O(n^), where oj < 3 is the exponent for matrix multiplication. But we only need to compute 
it once and save it. The rest of Step 1 takes 0{n 2 ) operations. 

Remark 3.13. It is possible to give bound on the running time that involve the diameter or the 
maximum effective resistance instead of algebraic connectivity, using the techniques appeared in 
[33]. [25] . We omit this discussion here. 
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3.5. Efficient bijective proof of the matrix-tree theorem. By Lemma [2T31 the order of Jac(G) 
is equal to the value of the determinant of the reduced Laplacian of G. In order to have a bijective 
proof of the matrix-tree theorem, one needs to find an explicit bijection between the elements of 
the group Jac(G) and the spanning trees of G. It is rather clear that such a bijection cannot be 
fully canonical, as that would imply a particular spanning tree is distinguished and corresponds to 
the identity! Therefore, one needs to make some choices to write down a bijection. 

By definition Jac(G) = Div°(G)/Prin(G). So the Jacobian of G is the set of equivalence classes 
of divisors of degree zero. We know by Proposition 12.51 once we fix a vertex q, there is a unique 
representative in each class that is g-reduced. 

Now, one can ask again whether there is a nice and explicit bijection between these g-reduced 
divisors and spanning trees of G. It turns out to have such a bijection, one needs to make a few 
more choices, for example choose an ordering on the edges of the graph. The problem of giving 
explicit bijection between reduced divisors (G-parking functions) and spanning trees of graph is 
studied extensively in the literature (see, e.g., [121 OS CE1 E] ) • For our work, we use the bijection 
given by Cori and Le Borgne in [18]. In fact, a variation of their bijection algorithm can be viewed 
as the Dhar algorithm with some extra features. The running time of the algorithm is then roughly 
the same as the running time of the Dhar's algorithm. We omit further details in this paper. 

Remark 3.14. Some of the bijections in the literature are given between the g-critical configurations 
(or sandpile predictions) and spanning trees. For a fixed vertex q, (/-critical configurations provide 
another set of unique representatives for equivalence classes of divisors (see, e.g., [3 [10]). There is 
a rather simple relationship between reduced and critical divisors; they add up to deg(v) — 1 for 
any v ^ q, and their value on q is forced by the degree of the divisor (see [2U [HI [HI 131] 15]). 

We have shown in §3.21 that this unique representative can be computed efficiently. Therefore we 
have the following "constructive and efficient" version of Kirchhoff 's matrix-tree theorem. 

Theorem 3.15. Size of the group Jac(G) is equal to the number of spanning trees ofG. Therefore 
the value of the determinant of the reduced Laplacian of G gives the number of spanning trees ofG. 
Moreover, there exists an efficiently computable bijection between elements of Jac(G) and spanning 
trees of of G. 

Appendix A. Some applications 

A.l. Uniform spanning trees. First application of our results is to give a new, completely al- 
gebraic method for picking uniform spanning trees of a given graph G. The running time of this 
algorithm is bounded deterministically. Proving that the output of the algorithm is indeed (uni- 
formly) random comes for free. The new algorithm is taking full advantage of the Jacobian of the 
graph which, as we observed in §2.31 is a natural group attached to the set of spanning trees of the 
graph. 

Here is the algorithm for picking a (uniformly) random spanning tree. Proof of correctness is 
trivial. 

Algorithm. (Picking a uniform spanning tree) 
Input: A graph G. 

Output: A uniform spanning tree of G. 

(1) Fix an arbitrary vertex q £ V(G), and an arbitrary ordering on E(G). 

(2) Compute the invariant factors {n\, . . . ,n s } and the corresponding generators {gi, . . . ,g s } 
for Jac(G) using a Smith normal form algorithm on Q (gj's are presented by divisors of 
degree zero). 
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(3) For 1 < i < s, pick a random integer < dj < n— 1, and compute the divisor D = J2t=i a iSs- 

(4) Use Algorithm 13.21 to find the unique q-reduced divisor in the equivalence class of D. 

(5) Use the Cori and Le Borgne algorithm in [18] (or its variation) to find the corresponding 
spanning tree. 

To our knowledge, the fastest Smith normal form algorithm given to date is given in [25], which 
runs in (n 2,697263 log ||Q||) 1+o( - 1 ^, where \\Q\\, for our application, means the maximal degree of a 
vertex. Note that for repeated sampling, one can perform steps (1) and (2) only once, and save 

invariant factors and generators. Step (3) takes O(-^p) chip-firing moves. Step 5 can be done in 
less than 0(n 2 ). 

With this algorithm, unlike previous methods, it is very easy to sample multiple spanning trees 
with certain joint distribution. For example, very few random bits are required to generate the 
pairwise spanning trees. We omit the details here. 

The bound on the running time of our algorithm does not beat the current best known running 
time 0(11^) of [17] . For repeated generation of spanning trees, it is not clear how our running 
time O^^/nmjX'i) compares with 0{n u ). We are optimistic that there is room for improvement in 
Algorithm 13.21 and/ or its running time analysis. Once again we emphasize that the techniques used 
for bounding the running time of Algorithm 13.21 is blind to the algorithm itself, and in practice the 
algorithm seems to run much faster than what is predicted by the upper bound. 

A. 2. The group law on reduced divisors or critical configurations. Fix a vertex q. Jac(G) 
induces a group structure on the set of g-reduced divisors (G-parking functions) or g-critical divisors 
(or sandpile predictions) of G. The latter is called critical group (or sandpile group) of G. The 
group law requires to first add the two given divisors as elements of Div(G), and then find the 
unique q-reduced or g-critical divisor equivalent to the addition result. Clearly our algorithm can 
be used to perform the group law. A different approach for performing the group operation is given 
in [3l]. Our algorithm seem to be faster, and our proof is significantly simpler and more algebraic. 

A. 3. Algorithm for the chip-firing game of Baker and Norine. The following Proposition 
is implicit in the proof of Theorem 3.3 in [5]. Proof of part (a) is, of course, trivial. 

Proposition A.l. Assume the chip-firing game of Baker and Norine starts with a configuration 
D 6 Div(G), and let D' ~ D be the corresponding q-reduced divisor. 

(a) If D'(q) > then there exists a wining strategy, and D' is a wining configuration. 

(b) If D'{q) < then no winning strategy exists. 

As a corollary, if one computes the (/-reduced divisor associated to a given configuration, then 
one can efficiently decide whether there exists a wining strategy or not. Moreover, it is a simple 
linear algebra exercise to find the moves for the winning strategy, once the initial configuration and 
the wining configuration is known. 

A. 4. Checking whether r(D) > 0. Baker and Norine prove a Riemann-Roch theorem for finite 
graphs in [5j. To any divisor D £ Div(G) they associate an integer, called rank, r{D) > —1. We 
refer to their paper for details. The following lemma follows from the proof of Theorem 3.3 in [5]. 

Lemma A. 2. r(D) > if and only if the unique q-reduced divisor equivalent to D has non-negative 
number on q. 

In fact, we can use this lemma and our algorithm to check whether r(D) > c for any "constant" c. 
Question of finding reduced divisors was first posed by Henrik Lenstra to Baker and Norine through 
personal communication. This connection with Riemann-Roch theory was the original motivation 
of the author to work on this problem. 
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